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PERCIVAL’S CONJECTURE FOR THE BUNIMOVICH 
MUSHROOM BILLIARD 

SEAN P GOMES 


Abstract. The Laplace-Beltrami eigenfunctions on a compact Riemannian man¬ 
ifold M whose geodesic billiard flow has mixed character have been conjectured 
by Percival to split into two complementary families, with all semiclassical mass 
supported in the completely integrable and ergodic regions of phase space respec¬ 
tively. In this paper, we consider the Dirichlet Laplacian on a family of mushroom 
billiards Mt parametrised by the length t G (0,2] of their rectangular part. We 
prove that Mt satisfies Percival’s conjecture for almost all t G (0,2], hence pro¬ 
viding the first example of a billiard known to satisfy Percival’s conjecture. 
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1. Introduction 


The Bohr correspondence principle informally asserts that the evolution of a quan¬ 
tum mechanical system coincides with the evolution predicted by classical mechanics 
in the large scale limit. 

One of the settings in which aspects of this correspondence can be made rigorous 
is that of dynamical billiards, which we shall now outline. 

If (M, g) is a compact boundaryless Riemannian manifold, we define dynamical 
billiards on M to be the Hamiltonian flow (pt on the cotangent bundle T*M of the 
manifold given by Hamilton’s equations 


( 1 . 1 ) 


dH 


dH 


~ Cl,- > ~ 


di 


dxj 


for the Hamiltonian H{x,^) := |(x,^)|^_i where g 


-1 


is the dual metric tensor. 

Since the Hamiltonian is a constant of motion for the flow </>t, it is natural to 
restrict the domain of this flow to the cosphere bundle 


(1.2) S*M := {z = {x, 0 G T*M : \z\g-i = 1}. 

More generally, one can define billiards on compact Riemannian manifolds with 
piecewise smooth boundary in the sense of Chapter 6 of [H], see also m 

To be precise, we assume that we can smoothly embed M in a boundaryless man¬ 
ifold M of the same dimension and that there exist finitely many smooth functions 
fj G C°°{M) such that the following conditions are satisfied. 

(1) dfi\ j-i(o) 


(2) dfi,dfj are linearly independent on ^(0) n fj ^(0), 


(3) M = {x € M : fj{x) > 0 for all j}. 

We can then write 

dM = UjdMj := n M) 

and denote by 5 C dM the set of points that he in dMj for multiple j. 


We dehne the broken Hamiltonian flow pt on S*M locally by extending the bound¬ 
aryless Hamiltonian flow by reflection at non-tangential and non-singular boundary 
collisions. 

That is, if ptoi^) = (^^j'^-i-) with x G dM \ S and > 0 where is 

the outgoing unit normal covector, we extend pt to sufficiently small t > to by 
defining pt{z) = pt-toix,^-), where G S*M is the unique covector such that 
G T*dM and 7r((^_|_) = vr(,f_) where vr : Tq^M —)• T*dM is the canonical 
projection. We terminate all trajectories that meet dM in any other manner. 
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There are four subsets of phase space for this class of manifolds which 

present an obstruction to obtaining a globally defined broken Hamiltonian flow or 
to the application of tools from microlocal analysis. We enumerate these sets below. 

(1) Bi = {ze S*M : e 5} 

(2) B 2 = {z € S*M : (j)t{z) G dM for infinitely many t in a bounded interval} 

(3) B 3 = {z G S*M : (j)tlz) ^ dM for any f > 0 or (ptiz) ^ dM for any f < 0} 

(4) Bi = {z G S*M : (ptiz) meets dM tangentially for some t G M.} 

Removing these sets from our flow domain, we then obtain a globally defined bil¬ 
liard flow on 2? = S*M\ (Wj^^Bj). For manifolds without boundary, we simply take 
V = S*M. 


The canonical symplectic form A dx on T*M determines a family of measures 
/ic on each of the energy hypersurfaces 

(1.3) T,c = {z = (x,^) G T*M : \z\g-i = c} 

defined implicitly by 

fd^cdc= / f\d^/\dx\ 

for / G Cf{T*M). 

Upon normalisation of we then obtain the Liouville measure /jl on S*M, which 
allows us to study the ergodic properties of the billiard flow cpt. 

Remark 1.1. It is shown in Section 6.2 of [6] that the sets Bi,B 2 are of Liouville 
measure zero, and it is shown in m that the set B 4 is of Liouville measure zero 
for the class of manifolds considered. That the remaining set B 3 is Liouville null 
is usually taken as an assumption. In particular, it is clear that this assumption is 
satished by bounded domains in M"". 



The billiard M is said to be ergodic if for /i 2 .-almost all z G 2? and every 
/ii-measurable A <G S*M we have 


(1.5) 


lim 

T^oo 


\{t G [0,T] : M^) ^ A}\ 
T 




The most famous example of an ergodic billiard is the stadium in M^, the ergodicity 
of which was first studied by Bunimovich [3]. 

Another prototypical example of ergodic billiards is provided by surfaces of con¬ 
stant negative curvature, where ergodicity is a consequence of the hyperbolicity of 
the flow. A proof of ergodicity in this setting can be found in Hopf m- 

On the other hand, if the billiard flow on M is completely integrable, then individ¬ 
ual trajectories are constrained to n-dimensional Lagrangian submanifolds specified 
by the values of the n constants of motion, and certainly do not equidistribute. 
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The quantum mechanical analogue of the system (1.1) is the evolution of a wave 
function ■0 G Lp‘{M) according to the rescaled Schrodinger’s equation 


( 1 . 6 ) 


- A„0 = i 


.dip 


dt 


with boundary conditions to ensure self-adjointness of the Laplacian. We shall 
work with the most studied and technically easiest choice of Dirichlet boundary 
conditions. 

Since the boundary of M is Lipschitz, it follows that the Laplacian —Ag is self 
adjoint on when given the standard domain {M) r\H q{M). Standard spectral 
theory then shows that — A^ has purely discrete spectrum (counting multiplicity) 
{0 < El < E 2 <...} C M+. 

By choosing a corresponding orthonormal basis of eigenfunctions {uj)j£fq, we can 
thus separate variables and formally expand solutions to (1.6) as 


(1.7) 


u{x,t) = ajUj{x) ex.p{—iEjt). 
i=i 


From this equation, we see that the localisation properties of high energy solutions 
to (1.6) are encoded in the high energy eigenfunctions of the operator — A^. 


The phase space localisation of the high energy eigenfunctions of — A^ can be 
described using the calculus of semiclassical pseudodifferential operators, as defined 
in Chapters 4 and 14 of m- 

To each subsequence of (uj), we can associate at least one non-negative Radon 
measure /i on S*M which provides a notion of phase space concentration in the 
semiclassical limit. 

We say that the eigenfunction subsequence (uj^) has unique semiclassical measure 
fj, if 


( 1 . 8 ) 


^Im (a(a;, E^^^D)uj^ , uj^) 



a{x, ^) dfi 


for each semiclassical pseudodifferential operator with principal symbol a compactly 
supported supported away from the boundary of S*M. In Chapter 5 of HZ!, the 
existence and basic properties of semiclassical measures are established using the 
calculus of semiclassical pseudodifferential operators (see also El)- 


A billiard M is then defined to be quantum ergodic if there is a full density 
subsequence of eigenfunctions (un^) such that the the Liouville measure on S*M is 
the unique semiclassical measure associated to the sequence Un^- This statement 
can be interpreted as saying that the sequence of eigenfunctions equidistributes in 
phase space with the possible exception of a sparse subsequence. 

It is a celebrated result due to Gerard-Leichtnam [8] and Zelditch-Zworski [16] 
that compact Riemannian manifolds that have ergodic geodesic flow are quantum 
ergodic. This generalises the earlier results of Schnirelman |13) . Zelditch m and 
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Colin de Verdiere [5] in the boundaryless setting. 

In this paper we consider the family of mushroom billiards Mt = Rt'<J S C 
where Rt = [—rijri] x [—1,0] and S is the closed upper semidisk of radius r 2 > ri 
centred at the origin. We denote the area of Mt by A{t). 



Figure 1: The half-mushroom billiard, with a high energy eigenfunction that 
extends by odd symmetry to the mushroom billiard. This particular eigenfunction 
appears to live in the ergodic region of phase space. Image courtesy of Dr Barnett. 


This billiard, proposed by Bunimovich [1] is neither classically ergodic nor com¬ 
pletely integrable for t > 0 and is rather one of the simplest billiards that satisfies 
the following mixed dynamical assumptions. 

• M is a smooth Riemannian manifold with piecewise smooth boundary 

• The flow domain T> is the union of two invariant subsets, each of positive 
Liouville measure and one of which, U, has ergodic geodesic flow 

• The billiard flow is completely integrable on 'D\U. 


In the mushroom billiard, Ut consists of /ii-almost all trajectories that enter 
Rt U B{0,ri) before their first boundary collision. The trajectories that do not 
enter Rt U B{0, ri) before their first boundary collision lie entirely within the upper 
semi-annulus S\B{0, ri) and are just reflected trajectories of the disk billiard. The 
integrability of the geodesic flow on T>\Ut then follows from the integrability of the 
disk billiard. 
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In the case of such mixed systems, we do not yet have a satisfactory analogue to 
the quantum ergodicity theorem. It is a long-standing conjecture of Percival |12j 
that a full density subset of a complete system of eigenfunctions of the Laplace- 
Beltrami operator can be divided into two disjoint subsets, one corresponding to 
the ergodic region of phase space and the other corresponding to the completely 
integrable region. Moreover, the natural density of these subsets is conjectured 
to be in proportion to the Liouville measures of the corresponding flow-invariant 
subsets of V. 

Conjecture 1.2 (Percival’s Conjecture). For every compact Riemannian manifold 
M such thafD is the disjoint union of two invariant subsets U,'D\U, with U ergodic 
and V\U integrable, we can find two subsets A,BcN such that 

(1) AL) B has density 1 

(2) {uk)keA equidistributes in the ergodic region U 

(3) Each semiclassical measure associated to the subset B is supported in the 
completely integrable region T)\U 

(4) The density of A is equal to yLi{U). 


Numerical evidence due to Barnett-Betcke [2] has strongly supported this con¬ 
jecture for the mushroom billiard. In this paper, we prove Conjecture |1.2| is indeed 
true for the mushroom billiard, at least for almost all t G (0,2]. Essential in our 
work is the following result due to Galkowski [7|. 

Theorem 1.3. For any compact Riemannian manifold with boundary satisfying Q, 
there exists a full density subsequence of (uj), such that every associated semiclas¬ 
sical measure /r satisfies 

(1.9) t\u = ap.L\u 

for some constant a. 

Our strategy for this proof is motivated by that used by Hassell in constructing 
the first known example of a non-QUE ergodic billiard [TO]. 


We begin in Section 2 by using the Dirichlet eigenfunctions on the semicircle to 
construct a family (un,a^) of 0 {n~°°) quasimodes that are almost orthogonal and 
are microlocally supported in the completely integrable region S*Mt \ Ut- 

Using the well-known asympotics of the Bessel function zeroes, we obtain a lower 


bound (2.11) for the counting function of this quasimode family. 


In Section 3, the main result is Proposition 3.1 an abstract spectral theoretic 


result that allows us to approximate certain eigenfunctions by linear combinations 
of quasimodes of similar energy given that the numbers of each are comparable. 
This is the essential ingredient for passing from localisation properties about our 
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explicit family of quasimodes to localisation properties of a family of eigenfunctions 
with asymptotically equivalent counting function. 


In Section 4 we commence our study of the variation of eigenvalues as the stalk 
length t varies in (0,2], In order to simplify the nomenclature, we often interpret t 
as a time parameter. 

The Hadamard variational formula asserts that 


( 1 . 10 ) 


Eit) = - Pt{s){dnu{t){s))‘^ ds 


IdMt 


where pt{s) is the unit normal variation of the domain at a boundary point s. For 


normally expanding domains such as ours, (1.10) directly implies that individual 
eigenvalues are non-increasing in t. 

However, using an interior formulation of the Hadamard variational formula from 


Proposition 4.1, we can also quantify the variation of the eigenvalue Ej{t) by 
(1.11) E-^{t)Ej{t) = {Quj{t),uj{t)) 

for an appropriate pseudodifferential operator Q supported in the stalk Rt C Mt. 

Proposition |4.2| then establishes that for a full density subset of the eigenvalues, 
the quantity {Quj{t),Uj{t)) can be approximated up to an error of o{Ej) by cutting 
off Q sufficiently close to the boundary dMt. This result is shown by using analysis 
of the wave kernel to establish the key spectral projector estimates ( |4.9[) a nd (4.10). 
We can then use the equidistribution result of Galkowski’s Theorem 1.3 to asymp¬ 


totically control {Quj{t),Uj{t)) and hence provide us with an upper bound (4.17) 
on the speed of eigenvalue variation for almost all eigenvalues. 


Section 5 completes the argument in two parts. 

In the first of these parts, we define a set G C (0, 2] such that for t € G, we have 
a certain spectral non-concentration property on Mt. Precisely, we have that 

the number of eigenvalues lying in the union [Oj — c,a'j + c] 
can exceed n by at most a small proportion, for large n. 


( 1 . 12 ) 


Proposition 3.1 then implies that these eigenfunctions are asymptotically well- 
approximated by linear combinations of the previously constructed family of quasi¬ 
modes {vn) which are microlocally supported in the completely integrable region 
S*Mt \Ut of phase space. 


In fact, the explicit computation (2.11) of the counting function of these quasi¬ 


modes leads to a proof that the corresponding family of eigenfunctions must fill up 
phase space. We show this is Theorem 5.3 


Consequently, we show in Proposition 5.4 that a full density subset of the com¬ 
plementary family of eigenfunctions must have all semiclassical mass in the ergodic 
region Ut- From Theorem|1.3[ this family must then equidistribute in Ut as required. 


The final part of the paper establishes via contradiction that (0, 2] is Lebesgue- 
null. As in m we can choose the eigenvalue branches Ej{t) to be in increasing order 
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and piecewise smooth in t. The crucial ingredient here is then the asymptotic bound 


(4.17) on the speed of eigenvalue variation. 


If Q is not of full measure, we can construct a small interval X = [ti,t 2 ] in which 
the average number of eigenvalues Ej (t) lingering near quasi-eigenvalues af exceeds 


d = d{ti) = 1 — fiiiUti) by using the negation of (1.12). 
Now Weyl’s law 

X^\Mt\ 


(1.13) 


Nt{X^ 


dvr 


implies that the decrease of eigenvalues over X is asymptotically given by 
(1.14) E,{ti) - E,{t2) ~ A7rj{A{ti)-^ - ^(ta)”^) 

in I as j ^ oo. 


We can use (1.14) together with the fact that the small windows about quasi¬ 


eigenvalues are comparatively sparse in the interval [Ej{t 2 ), Ej{ti)] to show that the 
upper bound (4.17) on eigenvalue speed provides a lower bound of (1 — d) on the 
time they must spend travelling outside of quasi-eigenvalue windows. 

This implies that the average proportion of time spent by large eigenvalues lin¬ 
gering near quasi-eigenvalues for t G X cannot exceed d, and consequently that the 
proportion of lingering eigenvalues cannot exceed d. This contradiction concludes 
the proof. 

I would like to thank my doctoral supervisor Professor Hassell for suggesting this 
problem and for onr many fruitful discussions regarding it. 


2. Quasimodes 

In polar coordinates, the Dirichlet eigenfunctions for the semidisk are given by 

( 2 . 1 ) Un,k ■= sm{n0)Jn{an,kr/r2) 

where an,k is the k-th positive zero of the n-th order Bessel function J„. 


Proposition 2.1. If we define 
(2.2) Vn,k 

where 


Xir)un,k 

IIX'^n,fc||L2 


(2.3) 


X{r) 



for r < ri 

for r > (ri -|- e)\/l — > n. 


then the family 

/ \ r / 9 / 9\ ^^2 T 

(2-4) {{Vn,k, Ot^,k/r2) ■ an,k < 


forms an 0{n °°) family of quasimodes, with all semiclassical mass contained in the 
completely integrahle region S*Mt \ Ut of the billiard. 
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Moreover, these quasimodes are almost orthogonal, in the sense that 
(2.5) \{vn,k,Vm,i)\ = 


Proof. The restriction on k in our family implies that the error incurred in cutting 
off only depends on the values of the Bessel function Jn{x) for x G [0, ny/l — e^]. 
Then from [T], we have the estimates 


( 2 . 6 ) 

and 

(2.7) 


Jn(nx)\ < - , — for X < 1 


xv^{i + vi^r 


for X < 1 


for bounding the Bessel function near 0. 

Together these estimates imply that the error incurred by cutting off is 0{n~°°). 
Furthermore, as the u^^k are pairwise orthogonal, these bounds also show that the 
Vn,k are almost orthogonal in the sense claimed. 


Now for any smooth compactly supported symbol a spatially supported in Rt U 
B(0, ri), the disjointness of supports from our family of quasimodes implies that 

(2.8) (a(x, ir 2 /an,k)D)vn,k,Vn,k) = 0{n~°°). 

In particular, we have that any semiclassical measure /x associated to these quasi¬ 
modes cannot have mass in the region {(x,^) G S*Mt : x £ RfU B(0,ri)}. 

Moreover, by the flow invariance of semiclassical measures (See Theorem 5.4 in 
m), this implies that any corresponding semiclassical measure cannot have mass 
in the ergodic region Ut because the pre-images under geodesic flow of 
{(x, ^) £ Vt ■ X £ RtD B{0, ri)} cover Ut- □ 

Proposition 2.2. ITe can index these quasimodes as (x„,a^) so that the quasi¬ 
eigenvalues are in increasing order, whilst having 

(2.9) (A + al)vn = 0{n-n = 0(«n~) 
and 

(2.10) \{vn,Vk)\ = k)~°°) = 0(min(Q;n, afc)““). 

Moreover, as e —>• 0, the counting function of these quasimodes has the following 
asymptotic bound. 


Proposition 2.3. 

. ,#{{n,k):an,k/r2<\an,k<^J 
(2.11 iimini-- 7 ^-^— 

^ ^ \ 


hUUt) 

hUBt) 



4?! 


where A{t) is the area of the mushroom Mt. 
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Proof. To simplify our calculations, we scale fiL so that = 2TrA{t). 

From an arbitrary point (r, 9) in the annulus, the trajectories that never enter 
the stalk have measure (27r — 4sin“^(ri/r)) out of the full measure 27r of the unit 
cosphere at that point. 

Hence 


p7T pr2 

^l(A) -= / / r(27r - 4sin“^(ri/r)) <ird6» 

Jo Jri 


fri 

— 7r^(r2 — rf) — 47r / r sin“^(ri/r) dr 


r»’2 


I ri 


= 7r^r| — 27rr^ — 1 — 27rr| sin ^{C 


where C = r^jTx. 
This implies that 


( 2 . 12 ) 


1 - 


fJ.L{Ut) \ A{t)\ ^ _ r. 
^J-L{'Ut)) 47r 


= ^ I- 




To estimate the left hand side of ( 2.11| ), we use the leading order uniform asymp¬ 
totics for Bessel function zeros found in [T]. 

As n —)> oo, we uniformly have 

(2.13) an,k = nz{n~‘^^^ak) + o{n) 

where z : (—oo,0] —)• [l,oo) is defined implicitly by 


(2.14) - 1 - sec-‘( 2 (C)) 

and the are the negative zeros of the Airy function, which have asymptotic 

(2.15) 0(VV3), 


We now write = r 2 / (ri -|- e). 


We count the left hand side of (2.11) by separating into two regimes based on the 
size of n/X. In each of these two regimes, a single one of the inequalities defining 
our family (2.4) implies the other. More precisely, we have 


\{{n,k) : an,k/r 2 < \ an,k < Cen]\ 

= \{{n,k) :n<r 2 X/Ce,an,k<Cen}\ 

+ \{{n,k) : r 2 X/Ce <n< r 2 X,an,klr 2 < A}| 

(2.16) =: iV^(A;e)+WB(A;e). 


For n, k sufficiently large, a sufficient condition for being in regime A of (2.16) is 
to have 


( 2 . 17 ) 


z{n ^^^ak) < Ce — e = Ce 
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and 


(2.18) n < raA/C,. 

Also, from the Airy function asymptotics we have 


(2.19) = A 

6 6n 




3/2 


_i, 

=-h 0 {n 

n 


Hence from the monotonicity of z, for all n, k sufficiently large with n < 
a sufficient condition for being in regime A of (2.16) is 


( 2 . 20 ) 


k < 


Ce — 1 — sec ^{Ce) — € 


-n. 


TT 


Noting that the contribution from small n and k is finite, we can conclude that 


lim inf — ^ > lim inf (—j 

A ^ 


A2 


n) 


C, - 1 - sec-\C,) - e 


n<r2XICe 


TT 


riiyCe -l-sec ACJ-e) 


2C2 


TT 


Similarly, for sufficiently large n,k, a sufficient condition for being in regime B of 
( |2.16 ) is to have 

(2.21) z{n~‘^^^ak) < \r 2 ln — e = A>e(A,n). 

and 


( 2 . 22 ) 


r 2 A/C'e < n < r 2 A. 


Hence, for all n, sufficiently large with r 2 XICe < n < r 2 A/(l + e), a sufficient 
condition for being in regime B of (2.16) is 


(2.23) 


k < - I - sec ^(A>e(A,n)) - 


TT 


Again throwing away a finite number of small pairs, we obtain 
Ns(A2) 


lim inf 

A A2 

1 


> lim inf E ( ny^D^iX, n)2 —1 —nsec ^{D^{X,n)) — enj 


Ce ^ - l + e 


ro A 

1 /■ife 

= lim inf —^ 

X vrA^ / r2^ 

Ce 


(^ty/Dt{X, ty — 1 — t sec ^{D^{X,t)) — et^ dt. 


Each of the three summands in the integrand has elementary primitive, so we can 
explicitly compute this quantity. 
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Noting that Ce,Ce —)■ C, we compute 

Na{X‘^) + Nb{X‘^) 


lim lim inf 
e-5>0 A 


A 2 


> 


+ 


r1\/ 6*2 — 1 rf 




271 

■y rr2\ 

ttX^ 




= 1 - 


r r tsec-(^) 

iriA V T^X^ Jr,\ t 

\/C2-l-^sin-i(C'-i)^ . 


dt 


71^2 


as required. 


□ 


3. Spectral Theory 

We next establish the following key spectral theoretic result. 

Proposition 3.1. Let LI be a Hilbert spaee. Suppose T G T(^) has a eomplete 
orthonormal system of eigenveetors {ui,Ei)i^^ with the sequence (Ei) non-negative 
and increasing without bound. Suppose further that we have a family of normalised 
quasimodes {vi, with 

(3.1) ||(T-ii;>,||<ei 
and 

(3.2) \{vi,vj)\<e 2 fori^j 

for some positive ei, 62 > 0. 

We write 

(3.3) V = Span{vi}^^i 
and 

n 

( 3 . 4 ) U = Span{uj : Ej G \^[E[ — c, E[ + c]}. 

i=\ 

We denote the orthogonal projection onto a subspace S CH by tts- 
If for some c > 0 and some 0 < e,S < 1/2 we have 


(3.5) 

and 


m = ff{j gN: Ej G[j[El- c, E'^ + c]} < n(l + e) 
i=l 


e 2 A 

(3.6) ^ + €2 < - 

n 

then at least n(l — ^/e) of the corresponding eigenveetors Ui satisfy 
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(3.7) \\ui - 7ry(iij)|| < + 2(5^/^. 

Proof. The idea behind the proof of the estimate ( |3.7| ) consists of several successive 
approximations. 

We first show that the projections Trjjivi) are almost orthogonal and can be trans¬ 
formed into an orthonormal basis of their span by a matrix A that is ap¬ 

proximately the identity. 

We then show that excluding some exceptional eigenvectors, the remaining eigen¬ 
vectors are necessarily rather close to the space W. This implies that the non- 
exceptional eigenvectors can be well approximated by their projections, which leads 
us to conclude u ^ t^wW) = Bw = BAttij{v) ^ BAv for some matrix B. 


To begin, we reindex the eigenpairs {ui,Ei) so that Ej G — c,E[ -|- c] 

precisely for j = 1, 2 ,..., m. 

The assumptions (3.1) and (3.2) then imply 


and 


\\{T - E'f)'^{vi,Uj)ujf 
jeN 

OO 

^ ^ \E, - E[Wvuu,)\^ 

j=m-\-l 

OO 

^ ^ \{Vi,Uj)\^ 

j=m-\-l 




< 


e 


2 

1 


< 


e 


2 

1 


< 


e 

c 


2 

1 

2 


> 



\{-Ku{Vi),'Ku{Vj))\ < \{Vi,Vj)\ + \{-Kur{Vi),-Kjjr{Vj))\ 

< £2 + Y^(l - ||7rr/(Ui)||2)(l - \\Tlu{VjW) 

< ^2 + 4 

& 


for i 7 ^ j. 


Together with (3.6) we obtain 


(3.8) \\Mv^)f>l 

and 


5 

n 


(3.9) 


\{7ru{vi),7ru{vj))\ < 


5 

n 


for i ^ j. 

The matrix M with entries Mij = {'n'uivi),7ru{vj)) satisfies 
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(3.10) 


||Af - I\\hs =■ II-E'IIhs < 6 < 1/2. 


Note that if the collection {Trjjivi)} were linearly dependent, then the matrix M 
would be singular. The estimate ( 3.10| ) precludes this possibility, because we can 
invert M = I — [I — M) as a Neumann series. In particular, this implies that m > n. 


We now write W = Span{7r[/(i;i)}/^^ and suppose that is an orthonormal 

basis for W which can be given by the transformation w = Attu{v), where A is an 
n X n real matrix that acts on the Hilbert space via matrix multiplication. 
Expanding out the matrix equation {wi,Wj) = 6 ij we obtain 


(3.11) AMA* = EE aikajl{'^uivk),7ru{vi)) = I 

k=l1=1 

which has a solution 

(3.12) A = M-^/2 ^ + ^ 


k=0 


k=0 


From (3.10), we then deduce 

OO 

(3.13) ||H - /||h5 < E < 2<5- 


k=l 


In the case m = n, that is when W = U, we can find a unitary matrix B with 
Bw = u. 

We now assume m > n, recalling that the assumptions of the proposition imply 
that this excess is small as a proportion of n. 

We have 

m m m n 

Y, \\Tiw^{ui)f = m-Y \Mui)f = m-YY \{ui, = m — n < ne 

1=1 1=1 i=l j=l 


which implies that 

(3.14) ■ hw^iui)f > V^} < ny/e 


and consequently 

(3.15) #{i : ||7rvi/(rij)|p > 1 — \/e} >m — ny/e > n(l — ^/e). 


We again re-index the eigenpairs for convenience, so that the first n' 
eigenvectors Ui satisfy the estimate in (3.15). 

In this case, we define the n' x n matrix B to have entries 


Ki -^^)1 


(3.16) Bij = {TTw{ui),Wj) 
and the vector u € by 

(3.17) {u.)ti. 
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We then have 

(3.18) tvl(w) = Bw 

and the i-th row Bi of B has norm trivially bounded by 1. 
This leaves us with 


(3.19) u = {u — 'Kwiu)) + BAt:u{v) = (m — 7rvu(u)) + BAv + BA{'ku{v) — v). 


which implies 


||(u - BAv)i\\'H 


< 

11^2 

:) + P^fi £2 ku(i’) - vWh 

< 

+ -{ 26 ) ■ 

■ Vd 

< 

,1/4 + 253/2. 



This estimate shows that each Ui has distance less than + 25^/^ to some 
element of V. 

Consequently 

(3.20) \\ui — 7ry(ui)|| < + 25^/^ 


as required. 


□ 


Our strategy to prove the main theorem is to control the number of eigenval¬ 
ues in most clusters formed by hnite unions of overlapping intervals of the form 
[a? — c, a? -|-c], and then repeatedly employ Proposition 3.1 to establish the exis¬ 
tence of a large density subsequence of these eigenfunctions that localises in the 
semidisk. 


4. Results on Eigenvalue Flow 


Central to the argument is the analysis of how eigenvalues flow as we vary t. 
Weyl’s law provides us with the asymptotic 


(4.1) 




AMt 

A-k 


where Nt is the counting function of the Dirichlet eigenvalues on Mt. 

To obtain a more precise statement about the change of individual eigenvalues, 
we employ an interior version of the Hadamard variational formula and Theorem 

oi 


In order to make use of Theorem 1.3, we choose (t>{y) G C“(M) non-negative, 


supported near y = —1/2, and with integral 1, and we define the family of metrics 
(4.2) gt = dx"^ + + l)(j)fdy‘^ 


on Ml. 

This metric induces a natural isometry It : {M\,gt) —)■ {Mt,gi). 

If we dehne Rt = {1 + {t — 1)/))“^/^, then we have the following result from 
Proposition 7 of the appendix of [TO] . 
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Proposition 4.1. Let u{t) be an -normalised real Dirichlet eigenfunction of A 
on Mt with corresponding eigenvalue Eft). We then have 


1 


Eft) = --{Quft),u{t)) 


(4.3) 

where the operator Q is given by 

(4.4) Q = -Adycftdy + [dy, [dy, (f>t]] = (ft - + (ftdl) 

on Mt- 

Here, (ft ■ Mt —>■ M is given by: 

(4.5) <ft = {(fRj)oI-\ 

We now cut Q off away from the verti cal s ides of the stalk so that we can use the 


interior equidistribution result Theorem 1.3 to control the quantity E, ^E^. 


We do this by defining 
(4.6) 

where X5 £ satisfies 


(4.7) 


x<5(a;) = 


Q& xsQ 


0 for X G [—ri, —ri + <5] U [ri — 6, n] 


1 for X G [—ri + 26, ri — 26\. 

Proposition 4.2. For any e > 0 and any t G (0, 2], there exists (5 > 0 such that 
(4.8) \Enl HQs - Q)Un^{t),Un^it))\ <€ 

for all k, where (n^) is a t-dependent subsequence of the positive integers with lower 
density bounded below by I — e. 

Proof. First we show that it suffices for each t to establish the spectral projector 
estimates 


(4.9) 
and 

(4.10) 


lhll[A,A+l)(V-^)llL2(Mt)^i“(Mt) — 0(A^/^) 


Here rjt = r] o Iff ^ where rj : Mi —>■ M is an fixed smooth cutoff function sup¬ 
ported and equal to 1 in a neighbourhood of dMi n spt{(f) such that r] vanishes in 
a neigbourhood of the semidisk. 

Applying r/tVlj;)^ ,)^+i)(-v/A) to X^a e[A A+i) using the estimate (4.10) then 

yields 


r]t{x) aj^Uj{x) 

Aj g[A,A+1) 



E 


ajUj 


Xjd [A,A+1) 


L2 

1/2 


E 


Aj G [A, A+l) 
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for each x G Mt- 

Setting Qj = Vuj{x) then yields the estimate 

(4.11) Y. |Vn,(x)|2 < 

Aj G [A,A+1) 

Similarly, we obtain 

(4.12) r}t{xf Y ^ 

Aj G[A,A+1) 


The estimates (4.11) and (4.12) then allow us to control each term of (4.8) in an 
average sense. 

For example, as only the the horizontal component 1 — of the cutoff function 
in (4.8) is J-dependent, we can integrate by parts in the second order term in (4.8) 
without loss. 

Then, by writing rjs to denote the cutoff function in the new second order term 
and choosing (5 > 0 sufficiently small so that rjs = rjijs, the contribution of these 
terms to (4.8) is controlled by 



E-^ 

E 

^7 



Ej<E 




Ei/2_ 

1 

r\j 

E-^ 

E 




k=l 




£l/2_ 

1 

< 

E-^ 

E 

k- 



k=l 





< 

CE- 


E 



\ 

k=l 

< 

Cs 




'M 


E-^ 


IM 


r]s{x)\'Vuj{x)\'^ dx 


Y ^ ^ [ V5{x)v{x) Y I 

k=l Aje[fe,fc+1) 


I Vuj(a:)p dx 



r]s{x) dx 


for sufficiently small 5 > 0, where —?• 0 as <5 —?■ 0. 

Together with analogous estimates for lower order terms, we obtain the estimate 


(4.13) 


1 


n 


YEYiiQs-Qhj,uj))\<Cs 




where —>• 0 as J —)• 0. 

By taking 6 sufficiently small that Cs < e^, we ensure that the collection of j with 
E~^\{{Qs — Q)uj,Uj))\ > e has upper density at most e. 


The estimate (4.9) follows from Proposition 8.1 in [9|. Note that the finite prop¬ 
agation speed of the operator cos(t\/—lA), the post-composition with a cutoff near 
a flat boundary, and the small-time nature of the argument together imply that Mt 










18 


SEAN P GOMES 


can be treated as the half-plane, which certainly satisfies the geometric assumptions 
of the cited result. 

By inserting the gradient operator in the dual estimate, it remains to control 
l|l[A,A+i)(\/^)V??f||Li(M0^L2(M0 in order to give us ( |4.10D . 

The argument in the proof of Proposition 8.1 in [9] allows us to replace the spectral 
projector by a smooth spectral projector A) where p'xis) = p{s — A) -|- p{—s — 

and p G 5(M) has non-negative Fourier transform supported in [e/2, e] for some suf¬ 
ficiently small e. 

So it suffices to estimate the norm of the operator 

(4.14) pTi^/^W = - [ cos(f\/AA)(e 

7T ./TO 


—itX^ 


X(t) -h e**^x(-f))Vr/(x) dt. 


This integral is supported close to t = 0, and hence by finite propagation speed, 
the kernel of the wave equation solution operator cos{t\/—A) on M is identical to 
that of the half-plane. 

Moreover, the kernel of the wave equation solution operator on the half-plane can 
be obtained from the free space wave kernel by the reflection principle, and their 
—)• norms are identical. 

This implies that it suffices to prove the estimate with the kernel for cos{t\/—A) 
replaced by the free space wave kernel. 

So the kernel to be estimated is 


Jr Jr2 Jr2 

= ^x{Kx{x,y))r]{x) 


where Kx,a\,'ilj are as in Lemma 5.13 from which we make use of in our 
penultimate line. 


Duality then completes the proof of (4.10) and the proposition. 


□ 


Proposition 4.1 allows us to use Theorem |1.3| and Proposition |4.2| to control the 
flow speed of a full density subsequence of the eigenfunctions for any fixed t. 


Proposition 4.3. For each t G (0,2] there exists a full density subsequence (n^) of 
the positive integers such that 

(4,15) liminfE„-‘{t)E„.(t) > 

where d(t) denotes the proportion of the phase space volume that is in the completely 
integrable region S*Mt \ Ut- 


Proof. From Proposition [L^ Proposition |4.2| and Theorem 1.3 for all e > 0 we may 
choose a d > 0 and a subsequence of eigenfunctions with lower density bounded 
below by 1 — e such that we have the estimate (|4.8[) and such that we have a unique 
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semiclassical measure jjL with = a^L\ut some non-negative constant a. We 
immediately have 


(4.16) 


a = 




< 


1 


1 — d(t) 1 — d(t) 
The definition of semiclassical measures then implies 

1 / 


liminf E (t)En^ > 
k-^-oo 


> -- 


> - 


2 Js*M 

1 


a{Qs)dn - € 


IS*M 


401^2 dn-e 


> - 


'i- — d Js*M 

1 


(t>t^2 - e 

■ 2 w 


> - 


Tr{l - d{t))A{t) Jm Jo 
A{ t) 


(j)t (x) sin^ {6)d6 dx — € 


— e. 


We can then apply Lemma |A.1| to obtain a full density subsequence of eigenfunc¬ 
tions satisfying the estimate (4.17). □ 

Moreover, we can strengthen the above to an almost-uniform result. 

Proposition 4.4. For any e > 0, there exists a full-density subsequence (n^) of 
positive integers and a family of sets Bnj. C (0,2] with m{BnA —>• 0 such that 


(4.17) 
for each t G (0, 2] \ B. 


Enlit)EnAi)> A{t){lAd(t)) ^ 


nk ■ 


Proof. For each <5 > 0 we define the subset 5^ C N as the collection of n G N such 
that 


(4.18) 


m{{t G (0,2] : E^^{t)Enit) < - 


Ait) 


A{t){l-dit)) 


- e}) > 6 . 


If every Ss were of zero density, we could write 5^ := N \ Ss and use Lemma A.l 
to assemble a full-density set satisfying the claims of the proposition. 

Now suppose that Ss has positive upper density for some <5 > 0. 

Since for every n G Ss, the sets B^ = {t : E~^it)Enit) < — ~ have 

measure bounded below and are subsets of a set with finite measure, there must 
exist a further positive density subset Ss C Ss such that 

Pi BnA^- 

nGSs 


(4.19) 
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This can be seen for instance by applying the bounded convergence theorem to the 
function 


(4.20) 


1 


n 




i=i 


The existence of a t in this intersection contradicts Proposition |4.3| and hence com¬ 
pletes the proof. □ 


The almost-uniform result in Proposition 4.4 can for our purposes be treated as 
a uniform bound on speed for large Ej, in light of the following weaker bound for t 


in the sets Bj of diminishing measure for which (4.17) does not hold. 


Proposition 4.5. There exists a positive constant C such that for every t G (0,2] 
and every j, we have 


(4.21) 


= -\{Quj,Uj) > -CEj{t) 


Proof. Integration by parts and Cauchy-Schwartz on the left-hand side provides us 
with a lower bound of 


-C 


M 


\{dyUj)uj\ -|- \{dyUj)uj\ + \uj\'^ dxdy 


> —C{{—Auj,Uj) + {—A.Uj,Uj)^^'^ + 1) 

> -CE, 

for a positive constant C that is uniform in time. 


□ 


Corollary 4.6. For any e > 0, there exists 5 > 0 and a full-density subsequence 
(ufc) such that we have 

(4.22) - ^ En, (t) dt < En, ih) (^2 - h) 

for any measurable set S C (0,2] with measure greater than 6. 


Proof. This follows from Proposition |4.4| and Proposition 4.5 by removing finitely 
many elements from the subsequence constructed in Proposition 4.4 □ 


We are now in a position to prove the main results of the paper. 


5. Main Results 

Definition 5.1. We call t G (0,2] good if for every e > 0, there exists some c > 0 
with 


(5.1) 


limsup+ ^1 + ^2 


n^oo ^ 

We denote the set of good times by Q. 
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I claim that Percival’s conjecture holds for the mushroom billiard Mt for all t G 
and moreover, that this set has full measure in ( 0 , 2 ], 

First we shall prove the claim for fixed t € G. 

For a given c > 0, we define c-clusters to be the connected components of 
Ui 6 N[af - c, + c]. 

We write A^semidisk(C))-^mushroom (C) to denote the number of quasi-eigenvalues 
and eigenvalues respectively contained in a given c-cluster C. 

The assumption then implies the following Proposition. 


Proposition 5.2. Suppose t and 0 < c < ‘llr\ are such that (5.1) holds, and all but 
finitely many c-clusters contain at least as many eigenvalues as quasi-eigenvalues. 
Then there exists a subset J of quasi-eigenvalues with lower density at least (1 — e) 
such that each quasi-eigenvalue in J is contained in a cluster C with 


(5.2) 


d^semidisk{C') ^ Tl^YiushroomiC) ^ (1 T c)Ng^^n.idiskif--') 


Proof. Index the c-clusters Ck in increasing order. 

Let 

S = {k : iVgemidisk(CA:) ^ -^mushroom(Cfc) ^ (1 T £)-^semidisk(l-'fc)} 

and 


^ {k . IVjnushroom(Cfc) ^ -^semidisk(Cfc)}' 

From the defining property ( |5.1[ ) of t ^ Q, we have: 

ON ,. -^mushroom(fc) 2 

(5.3) hm sup — - - - < 1 -h e 

n—>-00 / -‘^semidiskl^) 


The definition of S implies 


limsup 1 — 


'YLk<n -^semidisk(^)lF(^) 
Ylk<n -^semidisk(^) 


X)A:<n-^semidisk(^)(l If(^) ls(^))\ 2 

+ ^"-W- ITt -7 r\ - <l + e . 

Z^fc<n-^semidisk(fc) j 


The second term on the left-hand side is o(l) from the finiteness assumption. Hence 
we obtain that the upper density of N \ (5 U T) is bounded above by e. As F is 
finite, and consequently of density 0, we can conclude that the lower density of S is 
bounded below by 1 — e as required. □ 


Theorem 5.3 (Main Theorem). For each t € G, there exists Bt G N of density 
d{t) such that any semiclassical measure associated to the eigenfunctions {un)neBt 
is supported inside the completely integrable region. 


Proof. First, we fix e > 0 and choose c > 0 small enough so that the inequality (5.1) 
holds. 

Proposition |2.2| implies that we may choose the ei, 62 in applications of Proposition 
|3.1| to the increasing sequence of c-clusters to decay faster than any polynomial in 
the energy infima of these c-clusters. 
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By Weyl’s law, this ensures that for all but possibly hnitely many c-clusters, we 
have (3.6) with 5 decaying faster than any polynomial in energy. We remove the 
exceptional c-clusters, without any loss in density of our subset. 


In light of Proposition |5.2t we can then select a subset of the remaining c-clusters 
such that the included subset of quasi-eigenvalues has lower density exceeding 1 — e 
and such that ( |5.2[ ) holds for each cluster. 

We can now apply Proposition |3.1| on a cluster-by-cluster basis, with parameter 
(5 —)• 0 faster than any polynomial in energy. 

From the boundedness of pseudodifferential operators with compactly sup¬ 
ported symbols, this implies that for each e, we get a subsequence of eigenfunctions 
Ujf, such that any associated semiclassical measure /r satishes 

(5.4) fi{Vt \Ut)>l- 


Moreover, by comparison to Proposition 2.3 and Weyl’s law for the mushroom, 
we obtain a lower bound of d{t) — h{e) for the lower density of this eigenfunction 
subsequence with h{e) —)• 0 as e —)■ 0. So for each e > 0, we can find a c and 
a subsequence of {un) with density at least d{t) — h{e) which concentrates in the 
completely integrable region up to of its semiclassical mass. 

We now take a sequence Cj —)• 0 and denote the corresponding eigenvalue window 
widths by Cj. We write Bj^t to denote the corresponding concentrating eigenfunction 
subsequences. 

Lemma |A.1| then allows us to obtain a subsequence Bt of lower density at least 
d{t) such that any associated semiclassical measure /r satishes 


(5.5) ^,{Vt \ Ut) = 1. 

To bound the upper density of Bt, we choose a function Xe £ C“(M^ x M^) 
supported in the interior of M such that the following properties are satished. 

• 0 < X, < 1 


• Xe\vt\Ut - 0 

• !vt d^L > (1 - e)nL{Ut). 

Applying the local Weyl law (Lemma 4 from [16]) to the corresponding semiclas¬ 
sical pseudodifferential operator Xeix,hD), we obtain 

(5.6) 

^ {Xeix,E~^^‘^D)uj,Uj) > {l-e)nL{Ut) 

js[l,n]ni?t je[l,n]ni3f 


for sufficiently large n. 
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The localisation property (5.5) implies that the first summand is o(l) in n. Hence 
we have 

(5.7) ^ 


n 


{Xeix,E- ^^‘^D)uj,Uj) > (1 - 2e)^iL{Ut) 


j&[l,n]r\B^ 


for sufficiently large n 


From Theorem 


1.3 and the bound a < fJ-L{U) ^ that is immediate from semiclas- 


sical measures being probability measures, it follows that a full density subset of the 
remaining summands must be bounded above by 1 + e. 

This implies that 


(5.8) 


(1 + e) 


#{j <n:j€ Bf} 


n 


for sufficiently large n. 

Rearranging and passing to the limit n - 
required upper bound of 

#{j <n:j€ Bf} 


> (1 - 2e)iJLL{Ut) 


oo and then e —)• 0, we obtain the 


(5.9) 


lim sup 


n 


< 1 - miUt) = d{t). 


Hence Bt is a density d{t) sequence of eigenfunctions with semiclassical mass 
supported in the completely integrable region. □ 

Proposition 5.4. Let = N \ Bt. Then for each t ^ Q, a full density subsequence 
of {un)n&At equidistributes in Ut- 


Proof. From Theorem 5.3, the sequence of eigenfunctions {un)n&Bt has all semiclas¬ 
sical mass in the completely integrable region and Bt has natural density d{t). 

Applying the local Weyl law again with the function Xt from the proof of Theorem 
|5.3t we obtain 

(5.10) - 

n 


Y > l^L{Ut){l - e) 


je[l,npBf 

for sufficiently large n. 

Then, splitting the summation into the set 

(5.11) = {j E Bf : {xe{x,E~^^‘^D)uj,Uj) < 1 - ^/e} 

and its complement, the upper bound of 1 + e for a full density subset of the sum¬ 
mands in (5.10) then implies: 

(5.12) 


dn{A^,t) < hliUt) 


2 e 


e -h \/e 

using the notation dn from Lemma |A.1[ Passing to the limit n —?• oo, we obtain a 
subset of density exceeding 1 — 0 {^/e) of At with at least fJ.{Ut) > 1 — 0(\/e) for any 
corresponding semiclassical measure. 

An application of Lemma A.l then gives us a full density subsequence of At with 
all semiclassical mass in Ut. 
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Together with Theorem 1.3 this implies that we can find a full density subse¬ 
quence (rirtfe) of At such that every associated semiclassical measure is of the form 
'^Ut ■ fiL- □ 

We now show that the set (0, 2] \ ^ has Lebesgue measure 0. 

Proposition 5.5. If (0, 2] has positive Lebesgue measure, then there exists some 
e > 0 and some interval I = [^ 1 ,^ 2 ] C (0, 2] such that 


#{j G N : Ej{t) G ^2 ^ ^]| 


n 


dt > 1 + e 


(5.13) [ limsup 
l-^l Jx n—>-oo 

for all c > 0. Moreover, we can find such X with arbitrarily small length. 

Proof. By the monotone convergence property of measures, if m((0, 2] \ ^) >0 then 
there must exist e > 0 and a positive measure set S C (0, 2] on which we have 

(5.14) lim sup ( + >1 + 2 , 

n^oo \ ^ / 

for all t G S' and for all 0 < c < 2/r|. From the regularity of the Lebesgue measure, 
we can find an open set S C 17 C (0, 2] with m{U) < m{S) + 6 for an arbitrarily 
small 6 . We then have 

.,2 ^ ^,2 


(5.15) 
and 

(5.16) 


1^1 

1 


/ limsup 
'S 


/ #{j G N : E,it) G U(L,[af - c, of + c]} 




n 


dt > 1 - 1-26 


\u\s\ 


/ lim sup 
lu\S n^oo 


( #{j G N : Ej{t) G ur=iK^ - c, c]} 




n 


dt > 1. 


from our pointwise bounds on the integrands. By choosing 6 sufficiently small, we 
are thus guaranteed the estimate 

.,2 ^ „,2 

—t— t'\ ?- \ 

dt > 1 + e. 


1 


( #{j G N : Ej{t) G -c,af + c]} 


(5.17) — / limsup , 

1^1 Jjj n^oo \ 

Writing the open set 17 as a countable union of disjoint open intervals, the average 
of the integrand over one such interval must exceed 1 -|- e, as claimed. 

To complete the proof, we observe if we partition I into arbitrarily many intervals 
of equal length, at least one of them must also satisfy (5.13). □ 

To culminate the argument, we seek out a contradiction coming from the upper 
bound (4.17) on speed of eigenvalue variation and the lower bound (5.13) on the 
average proportion of eigenvalues lying in c-clusters. 

Proposition 5.6. For any e > 0, there exists c > 0 such that 


(5.18) 


lim sup — 

m^oo . T 


\{t € I: Ej € Ui [a? — c, of + c] 


< d(t\) + e 


for any sufficiently small interval X. 
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Proof. Note that we have the flow speed bound (4.17) for a full density subsequence 


of eigenvalues, so if we can establish the claimed inequality for each summand with 
a sufficiently large index that obeys the flow speed bound, density will allow us to 
draw the desired conclusion. 

We now suppose Ej is a large eigenvalue that lies in this full density subsequence. 

Writing X = — A{t 2 )~^) for brevity, Weyl’s law applied to the mushroom 

gives 


(5.19) 


Ej{ti) - Ej{t 2 ) > (dvrX - 2S)j 


for (5 > 0 and all sufficiently large j. 

Weyl’s law for the semidisk (recalling that we constructed the completely inte- 
grable region quasimodes from semidisk eigenfunctions) gives us an upper bound 
of 


(5.20) 


=|4 + <5 


for the number of quasi-eigenvalues in [Ej{t 2 ), Ej{ti)] and hence an upper bound of 


(5.21) 


2c 


TTr^X 


+ 5 j 


for the length of [Ej{t 2 ), Ej{ti)] that lies within Uj[a? — c,af + c]. 

Now suppose that Ej{t) spends proportion qj of t G X in — c, af + c]. 

From Proposition |4^ it follows that the qj are uniformly bounded above by some 


1-5. 


This means that we apply Corollary 4.6 to find a lower bound for the time taken 
by an eigenvalue Ej in our full-density subsequence to traverse the set 
[Ej{t 2 ), Ej{ti)] \u4af — c,af + c\. Heuristically, we can think of this as dividing 
the size of this set by an upper bound for the speed of the eigenvalue’s variation. 
Precisely, we have 


(5.22) {I - qj){t2 - ti) > 


jW(47r — 7rr|c) — j6{l + 2c) 


(5.23) 


(5.24) 

(5.25) 


Ejih 


> 


> 


X(47r — Trr^c) — <5(1 -|- 2c) 


4(ti) 


+ 5 


4-71 

A(ti) 




X(47r — 7rr|c) — <5(1 -|- 2c) 


4(ti) 

A(ti)(l-d(ti)) 


+ s 


XA(ti)^(l-d(ti)) 

A(ti) 


where the final two lines follow from Weyl’s law and passing to sufficiently small S 
and c respectively. 
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Additionally, since A{t) is a linear polynomial in t, we have 

A{t2) — A{ti) {t2 — ti)A{ti) 


(5.26) 

which implies that 


X = 


A{ti)A{t2) 


A{ti)A{t2) 


1-qj > 


liMn 

A{t2)^ 2 


qj < d{ti) + {d{ti) - 1) 
< d(ti) + e 


^ _ A 1 

A(t2) ) ^ 2 


for sufficiently small |X|, using the uniform continuity of A. 

Thus we have the required inequality for all sufficiently small intervals X and all 
sufficiently large j in a full density subsequence on which (4.17) holds. 

□ 


Proposition 5.7. For any e > 0 there exists c> 0, such that 


(5.27) 


7— / limsup 
Jt n^oo 


#{j £ N : Ej{t) £ - C, af + c]} 


n 


dt < 1 + e 


for all sufficiently small \X\. 


Proof. By the dominated convergence theorem, it suffices to show that we can find 
c such that 


J_ f #{j : Ej{t) £ - c, af + c]} ^ ^ ^ e 

\X\ Jt n 2 


(5.28) 

'■'lx n 

for sufficiently large n. This quantity is bounded above by 

1 |{t £ X : Ej £ Uj[Q;? — c, o? T c]}| 1 

n _ „ ” \I\ 


(5.29) 


j:Ej(t2)<al+c 


V 

n 

j-.Ej{t2)<al+c 


Qj- 


The sum is controlled by the previous proposition, giving us an upper bound of 

(5.30) — • max{j : Ej{t 2 ) < aj + c}{d{ti) + <5) 

for sufficiently large n. 

From Weyl’s law, we have 


(5.31) 


max{j : Ej{t 2 ) < aj + c} < {aj + 




for sufficiently large n. 

By taking 5, c, and I sufficiently small, we then obtain 
(5.32) i qj<' 

j-Ej{t2)<af+c 


n 


A{t2)d{t2)'^ 


dvr 
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Inverting the estimate (2.11) provides an upper bound of 1 + | for the first 
summand on the right-hand side for all sufficiently large n, thus completing the 
proof. □ 


Corollary 5.8. The set Q has full measure in (0,2]. 

Proof. This is an immediate consequence of Propositions 5.5 and 5.7. 

Appendix A. 


□ 


In this appendix, we prove the following abstract lemma that we have used several 
times to assemble full density subsequences along which a given function has limit 
0 . 

Lemma A.l. If there exists a function g : N —>■ M and a family of subsets Sj C N 
such that 


n^oo n 


(A.l) 
and 

(A.2) limsup g{n) < e'j 

nGSj^oo 

where Cj, e'- \ 0, then there exists a subset S CN such that 
(A.3) 




n 


and 

(A.4) 


lim g{n) = 0. 

nSS—>-00 


Proof. For ease of notation, we dehne 


(A.5) 


dn{A) = 


#{A: <n:kGA} 


n 


for A C N and n G N. 

We have g{n) < 2e'- for cofinitely many elements of Sj, and we denote these sets 
by Sj. Now let 

(A.6) Bj = {keN: g{k) > 2e' } C N \ 5'. 

Since each dn respects the partial ordering of set inclusion and is additive with 
respect to disjoint unions, we can construct a strictly increasing sequence {Nj)j^fq 
such that Ni = 1 and dn{Bj) < 1 — d + 2ej for all n > Nj. 

We define 

(A.7) 5= |JSjn[A(,-,oo). 

JSN 

If n G [Nj, Nj+i), then any k € [1, n] Pi i? must lie in Bi for some i < j and hence 
in Bj. 
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This implies that for n G [A^j, A^j+i) we have dn{B) < dn{Bj) < 1 — d + 2ej and 
consequently, that limsupdn(^) < 1 — d. 

71^00 

We now take S := N \ B, with the required density bound 


(A.8) 


liminf dn{S) > d. 

n^oo 


To complete the proof we observe that if n G [Nj, oo) H S, then n G N \ for 


each i < j, and hence g{n) < 2e'-. This establishes (A.4). 


□ 
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